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Abstract. A numerical study of the string tension and of the masses of the 
lowest-lying glueballs is performed in SU(iV) gauge theories for 2 < iV < 8 
in D=3+l and 2 < N < 6 in D=2+l. It is shown that for the string tension 
a smooth N — ► oo limit exists that depends only on the 't Hooft coupling 
A = g 2 N. An extrapolation of the masses of the lightest glueballs to A = oo 
using a power series in 1/A 2 shows that the leading correction to the infinite 
N value accounts for finite N effects for A at least as small as 3 and all the 
way down to A = 2 in many cases, fc-string tension ratios and possible issues 
connected with correlation functions at large A are also discussed. 



1 Introduction 

A better understanding of the large A limit of SU(A) gauge theories is an 
important prerequisite for advances in several branches of particle physics. It 
has been shown that accurate numerical results for observables in SU(oo) can 
be attained by extrapolating to A = oo values obtained at finite A ^ |2J E|. 
The extrapolation is based on the diagrammatic expectation that observables at 
finite A differ from the corresponding observables at A = oo by corrections that 
are described by a power series in 1/A 2 [I]. In this work we will present results 
for the string tension, /c-string tension ratios and glueball masses in 3+1 [3] and 
2+1 dimensions (Hj. 



2 The method 

We have studied numerically on the lattice SU(iV) gauge theories for 2 < N < 8 
in 3+1 dimensions and 2 < N < 6 in 2+1 dimensions. For any given N we have 
used the Wilson action 

s = f j E ( l - ^ Tr (tvw + ^(0)) , (i) 

l,fJL,V 
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Figure 1. In our scheme, a smeared link is obtained by adding to the original link (a) the 
parallel transported nearest neighbour links (b) weighted by a coefficient a and the next to 
nearest neighbour links parallel transported along all possible paths on the elementary cube 
(c) weighted by another coefficient /3. For a judicious choice of a and f3 it is possible to build 
operators with a 98% overlap on the physical states. 



where the U^y^i) are path ordered products of SU(iV) link variables along the el- 
ementary lattice squares (plaquettes) and (3 = 2N/(g 2 a 4 ~ D ), with g the coupling 
of the theory, D the dimensionality of the system and a the lattice spacing. The 
values of (3 used in our study have been chosen above the bulk phase transition 
point 0. 

Masses can be extracted from the exponential decay of correlation functions 
built with traces of path ordered products of functions of link variables. In order 
to have a significant overlap with physical states, these functions must be smooth 
on the characteristic distance scale of those states. Traditionally, iterated smear- 
ing and blocking |Sj are used to achieve this result. In our work in D=3+l we 
have used a combination of improved smearing (see Fig. ^) and blocking, while 
simple blocking has been used in D=2+l. For each state of interest, a set of oper- 
ators carrying the required quantum numbers has been built out of smeared and 
blocked variables and the corresponding correlation matrix has been computed. 
Masses have been obtained by a variational calculation performed on this matrix. 
More details on the numerical methods we have used can be found in [SJIHl- 



3 Correlators and large N 

There is a potential problem with the application of the method outlined in the 
previous section at large N. In fact, from general arguments, as N — > oo we 
expect factorisation of correlators, which implies that the connected part of the 
2-point function (from which we extract masses) is 0(1/N 2 ). At large enough N 
then this signal might become invisible. However, what determines the accuracy 
of the calculation is not the absolute magnitude of the signal, but the ratio signal 
over noise. Large N arguments 5 a show that this ratio is independent of N, and 
this expectation is confirmed by our data. Thus, despite the fact that the signal 
disappears as 1/iV 2 , we are able to compute correlators with an accuracy that is 
independent of N. 
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4 String tension 

The string tension a can be obtained from torelon correlation functions [2j. 
Diagrammatic arguments suggest that equivalent physics across SU(iV) gauge 
groups is obtained at fixed A = g 2 N. Since A depends on the scale at which 
the physics is probed, equivalent A are expected at the same physical scale, or 
equivalently the same scale is identified by a given A at any large enough value 
of N. Although this is not an ideal choice, on the lattice we can easily verify 
this statement on length scales of the order of a. At this scale lattice artifacts 
resulting from a naive definition of the lattice coupling can spoil the comparison. 
To limit the effects of those artifacts we have used the mean field improved 
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coupling f3j = [5/ < U^ u >, where < U^ v > is the average plaquette. Fig. El 
shows a plot of ay/a as a function of Xj = 2N 2 / f3i in D=3+l. This plot suggests 
that a smooth large ./V limit for \j exists at all values of a^fa. 

In D=2+l the lattice spacing a can be traded with the coupling g 2 , which now 
has dimensions of mass. Our data are displayed in Fig. El From diagrammatic 
arguments, we expect that up to 0(1/N 2 ) corrections, in the continuum limit 
yjojg 2 (x N. A fit gives 

n iQ7^o^ 0- 12 0(3) r9 s 
= 0-19755(35) ^— . (2) 
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Figure 4. Ratios of the tensions of strings of various A/"-alities over the tension of the funda- 
mental string. The scale on the y axis has been chosen in such a way that the values R — 1 and 
R — 2 always correspond respectively to Casimir scaling and to the sine formula. Those values 
are indicated by the horizontal lines in the plot. 



5 /c-strings 

Strings connecting sources in a representation of rank (or M-ality) k < N of 
SU(A^) with their antisources are known in literature as k-strings. Since strings 
associated to sources of different jV-ality can not mix, at fixed iV there are 
int[JV/2] stable /c-strings. The values of their tensions are determined by the 
underlying dynamics of colour confinement. 

/c-string tensions can be extracted from correlators of multi- wrapping torelons 
and powers of them that under the centre transformation e 2m l N acquire the phase 
e 2nlk / N . Our numerical calculations show that the ratios R of tensions of strings 
connecting sources of A/"-ality k and N — k over the tension of the fundamental 
string lay between the Casimir scaling (R = k(N — k)/(N — 1)) and the sine 
formula (R = sm(kir/N)/ sm(ir/N)) predictions in D=3+l [31 |H] (see also [TUj . 
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0++ 


0.8116(36) 


-0.090(28) 


3.28(8) 


2.1(1.1) 


++* 


1.227(9) 


-0.343(82) 


5.93(17) 


2.7(2.0) 


Q + + ** 


1.65(4) 


-2.2(7) 






0~ 


1.176(14) 


0.14(20) 






* 


1.535(28) 


-0.35(35) 






2++ 


1.359(12) 


-0.22(8) 


4.78(14) 


0.3(1.7) 


2 ++* 


1.822(62) 


-3.9(1.3) 






2— 


1.615(33) 


-0.10(42) 







Table 1. Fit parameters for the extrapolation to the N — oo limit of masses of glueballs 
in D—2+1 according to formula JIJJ and in D=3+l according to formula @. The quantum 
numbers of the states are indicated with the standard notation J PC and higher excitations in 
a given channel are denoted with a corresponding number of *. 

where a better agreement with the sine formula is claimed), while in D=2+l 
Casimir scaling seems to be favoured [§]. Although the corresponding numerical 
values are very close, those formulae hide very different scenarios: while the sine 
formula can be expanded in a series that contains only even powers of 1/N, the 
first term in the expansion of the Casimir prediction is proportional to 1/N, 
and this is hardly accommodated by diagrammatic arguments |111 I12j . Even 
if both formulae at small N are expected to receive sizeable corrections that 
undermine the numerical comparison, the question of the power dependence of 
the leading correction in a 1/N expansion still stands. Current numerical data 
are not accurate enough to provide a reliable answer to this problem. 

6 Glueball masses 

Glueball masses are extracted from correlators of products of smeared and 
blocked links along contractible paths. The transformation properties of the 
paths under rotation, parity and charge conjugation determine the quantum 
numbers of the states. We have fitted our data with the ansatz 

m N /g 2 N = do + djN 2 (3) 

in D=2+l and 

m N /ayfo = do + di/N 2 (4) 

in D=3+l. The results of our fits are reported in Table 1. In D=2+l, the ansatz 
works for all the states all the way down to N = 2. In D=3+l we draw similar 
conclusions for the ++ glueball, while a correction 0(1/N 2 ) turns out not to be 
enough to describe the SU(2) numerical data for the ++ * and the 2 ++ states. 
Note that while both in 2+1 and 3+1 dimensions we can extract accurate values 
for the adimensional masses in the N — > oo limit, the values of the coefficients of 
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the leading correction is determined unambiguously only in the D=2+l case. In 
D=3+l we can however estimate their order of magnitude. Our results support 
the conclusion that at finite N the deviation from the N — » oo value is a genuine 
0(1/N 2 ) correction. 

7 Conclusions 

Our investigation shows that a smooth large N limit exists for several physical 
observables and the physics at finite N is correctly described by the leading 
expected correction all the way down to at least N = 3. In particular, we have 
discussed our results for the masses of the lightest glueballs and for the string 
tension in both D=2+l and D=3+l. Unlike early works on this subject, in most 
cases we were able to provide accurate and reliable numerical estimates for the 
N = oo value and the leading correction. An interesting problem is whether for 
some observables the approach to the N = oo limit can be described by a 1/N 
correction. On the light of present results, this possibility can not be ruled out 
for ratios of /c-string tensions. 
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